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Abstract: The two-point gauge-invariant correlation function of gluonic field strengths, 

which is the main input in the stochastic vacuum model, is derived by using its relation 

to the Green functions of one- and two-gluon gluelumps. These Green functions are found 

by evaluating the path integrals of one and two gluons confined in the field of an infinitely 

heavy adjoint source. This allows one to obtain also the terms in the exponential asymp- 

totics of the two-point correlation function, which are produced by fluctuations of the 

gluons around their mean separation from the source. The pre-exponents are derived as 

well, demonstrating some novel dependences on the distance, 
c ■ 
• T—i - 

X 
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1. Introduction 

Despite a remarkable impact of the stochastic vacuum model 0, § on QCD phenomenology 
(see e.g. Q for a recent review), a field-theoretic establishment of this model is still pending. 
The main problem is an analytic calculation of the nonperturbative content of the two-point 
gauge-invariant correlation function of gluonic field strengths, which can be parametrized 
as § 

g 2 

— tr (F pu (x)$ xx <Fx p (x')$ x > x ) YM = (S^xSvp - 8 lxp 5 v \)D(x - x')+ 

+\ [dl ((z " x')x5 up - (x - x') p 5 vX ) + 8 X U {(x - x') p 5 pX - {x - x')a<W)] D^x-x'). (1.1) 

Here, & xx > = V exp [ig f®, dz p A^(z)T^ und j is a phase factor along the straight line, which 
goes through x' and x; r f " nd 's stand for the generators of the group SU(_/V C ) in the funda- 
mental representation; the average (. . .) Y m i s taken with respect to the Euclidean Yang- 
Mills action, \ f d 4 x(F^) 2 ; F* u = d p A a v - d„A» + gf ahc A\A%, a = 1, . . . , N 2 C - 1. 

Numerous lattice measurements of the correlation function (|l.l| ) were performing over 
the years (see jl], Q for the calculations with the cooling method and ref. HJ, where this 
correlation function was reconstructed from the lattice data on the static quark-antiquark 
potential). As for the papers attempting to analytically evaluate this quantity, some 
of those were dealing with the Abelian theories, where confinement was produced by 
monopoles |], f7|. At the same time, some progress was achieving in the non- Abelian 
case as well ||, |9|. In particular, a relation between the correlation length of the vacuum 
(that is the length at which the functions D and D\ experience an exponential fall-off) and 
the gluonic condensate was derived in || . Further analysis of the correlation lengths of D 
and D\ was done in Q . 1 Recently, a complete derivation of these functions was presented 



in [10]. Using initially the idea |ona relation between the correlation function (1J) and 



the Green functions of one and two (valence) gluons confined in the field of an infinitely 



1 \t was supporting the conjecture, tested on the lattice in pj, that these lengths might be different. 
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heavy adjoint source (the so-called one- and two-gluon gluelumps), the paper [10] further 
accommodated the Hamiltonian formalism for the analysis of these Green functions. This 
eventually yielded certain nonperturbative expressions for the functions D and D\. 

In the present paper, the Green functions (as well as the functions D and D{) will be 
calculated upon an evaluation of the corresponding integrals over trajectories of the gluons 
in the gluelumps. It turns out that, when the minimal-area-law Ansatz is applied to the 
Wilson loops of these (confined) gluons, further calculation of the path integrals becomes 
possible by making use of the Feynman variational method. By means of this method, we 
manage to find some novel terms in both the exponents and pre-exponents of the functions 
D and D±, which stem from the vibrations of the gluons around their mean separations 
from the heavy source. 

The paper is organized as follows. In the next section, we demonstrate our method 
at the example of a one-gluon gluelump and calculate in this way the function D±. In 
section 3, we proceed with the two-gluon gluelump and the function D. In section 4, a 
brief summary of the obtained results is presented. 

2. Green function of a one-gluon gluelump and the function D\ 

In the approximation when the gluon spin term is disregarded [10], the one-gluon gluelump 
Green function has the form 6^, U G(T, 0), where 

G(T, 0) = jds J Vx„ exp (- J* dr^A (W(C)) YM . 

x 4 (0)=0, x 4 (s)=T ^ ' 

Here, in the Wilson loop {W(C)) YM = tr (v exp (ig § c dx^A^T^^ , T a a dj 's stand for 

the SU(iV c )-generators in the adjoint representation, (T^) bc = —if abc ; the closed contour 
C is formed by the piece of the 4-th axis from the point (0, T) to the origin and the 
trajectory x^(t) of a gluon, which moves back from the origin to that point. Adopting for 
{W(C))ym the minimal-area-law Ansatz, we have 

oo 

G(T, 0) = N C (N* - 1) J ds J V Xi J Dxx 



x exp 



x 4 (0)=0, x 4 (s)=T 

dr I — + — + a ad j|x|x 4 



It will further be assumed that, in the directions transverse to X4, gluon, along most of its 
path, only slightly fluctuates around some mean distance, L, from the 4-th axis. 2 This 
means the following parametrization x(r) = L + l(r), where 1(0) = l(s) = 0, L = Ln. 
Then, up to the term linear in 1, |x| ~ L + nl. In this approximation, we have 



CO 



G(T, 0) ~ N c (Nl - l) e - a ^ LT J ds J Vx A j VVDfix 

2 4 (0)=0, Xi{s)=T 

2 This distance will be found numerically at the end of the calculation. 
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x exp 



dr - + ^ + £ + 



2,u 



where we have introduced an integration over the "einbein" field [i in order to make the 
potential term quadratic in 1. Next, the unit vector n may, without loss of generality, be 
chosen along the 1-st axis, n = (1, 0, 0). The integrations over liij) and h(r) then become 
those of a free particle. As for the integration over h(r), since 1, by definition, only weakly 
varies along the trajectory, one may replace l\ (r) by its value averaged along the trajectory, 
l\ — ► i Jq drlf. Altogether we have 



G(T,0) ~ iV c (iV c 2 - 1K 



°adj 



It I -/ 

X4(0)=0, x 4 (s)=T 



x exp 



dr 



| + f + l + u lK 



a 



adj 

2s 



dr^ 

M 



/ 2 



(2.1) 



Although the partition function of a harmonic oscillator, resulting now from the integration 
over h(r), is known exactly, its substitution would be of little use for a practical calculation. 
(Merely since the frequency of this oscillator is a functional of X4//X, and the subsequent 
integrals over X4 and /i can hardly be handled.) Rather, it is much more useful to apply 



Feynman variational method [11], that is fully legitimate, since l\ only slightly fluctuates 
around the origin. This method yields the leading quantum correction to the classical 
expression 



1 



st/(0) 



'4irs 
~ Jo dr 



+ U(x) 



(2.2) 

The correction originates 



for the integral J T>x{t) exp 

x(0)=x(s)=0 

just from the fluctuations of the trajectory around the origin and leads to the follow 



ing expression for the above-mentioned integral (see e.g. 12] for the details of derivation): 



/4-7TS 



exp — \ / — J_°° dxll(x)e « . The classical approximation ( |2.2| ) is recoverable from 

this formula in the limit s — > 0. 

Substituting for U the expression from eq. ( |2.1[ ), we then get for the path integral over 

l ^y- vfc ex p f-# Jo dT f) ■ This y ields 



G(T,0)=N c (N ( 



l) e -°-adj 



LT 



ds 



(4ns) 3 / 2 



I 



Vx 4 Vfix 



X4(0)=0, x 4 (s)=T 



x exp 



A.2 



dr ^ + Z + 



12 



(2.3) 



Further integration over \i replaces the sum ^ -\ — j^--^ by a single term o"adjv s/6|«4|. 
This is a world- line Lagrangian of a free scalar particle, which moves along the 4-th axis. 
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(The mass of this particle m = 0" a dj \/ s /6 i s > however, a function of s.) Such a Lagrangian 



7 adj 

is known to be equivalent to the one usually used for the world-line calculations, ^ + m 2 , 
so that the total world-line action reads Jq dr-^ + m 2 s. To restore the standard coefficient 
1/4 at the kinetic term, one can perform the rescaling r = 2r neW ) s = 2s new . The integral 
over X4 then becomes standard, and we obtain 



G(T,0) = N C (N ( 




To calculate the integral over s, let us change the integration variable to £ = (as 



1-1 



. The saddle-point approximation 



that yields for this integral a J °° d£ exp 
to the integral J °° d^e~^^ works well when the width of the corresponding Gaussian 



2 



distribution, which is 



(where £0 is the saddle point), is much smaller than 1. In 



our case, 



—7=^ = = (32/9) 1 / 3 (fj a djT 2 ) 2 / 3 , therefore the saddle-point approximation is 

legitimate when the distances under consideration are large in the sense of the inequality 
cr a djT 2 3> 1. Calculation of the integral is then straightforward, and we have 



G(T,0) =N c (N*-l) 



(WIT 

3 2 / 3 • 2 7 /3 . ^3/2 



exp 



The function D\(x) is related to this expression as [ 10 1 D\(x) 
the differentiation we obtain 



(9/32) 1 /3 ((JadjT 2 ) 2/3 

Performing 



2g 2 dG(x O) 
"W dx 2 ■ 



Dt(x) 



Cia, 



7T 



(32/9) 



1/3 °"adj 



x- 



1 L 

12V3 + 2 



°"adj 



1/3 



+ 



3(a adj x 2 ) 2 /3 



x exp 



-a adj L|a;| - (9/32) l/ *(a^ 



,2\2/3 



at (T a djX > 1, 



where Co 



AL 2 -l 



(2.4) 



2j y is the quadratic Casimir operator of the fundamental representation. 
Note that a gluon is confined only at large N c , where <7 a dj = 2<7f un d (while, at -/V c ~ 1, it 
is screened by other gluons). In the large- iV c limit, C2a s — * where A = g 2 N c is the 
so-called 't Hooft coupling, which stays finite in this limit. 

Equation ( |2.4D can further be compared to the expression for D\(x) found in ref. [10] 
by the Hamiltonian approach, which reads 



(2.5) 



where M ~ 1.5 GeV at <7 fund = 0.18 GeV 2 . One can see that the term (9/32) 1 / 3 (a a djX 2 ) 2 / 3 
in eq. ( [2.4D , although parametrically dominating the exponent in the limit cr a djX 2 3> 1, 
is numerically smaller than the leading term <7 a dj£|a;|. This is not surprising, since the 
term (9/32) 1 / 3 (<r a dj2; 2 ) 2 / 3 , originating from small fluctuations of the valence gluon around 
its mean separation from the heavy source, can only be a correction to the leading result. 

(9/32)1/3(^2)2/3 _ (9 CTf 2 und M) 1/3 



Indeed, equating a^L to Mq, we have 



T adj 



LlxJ 



2Mq 



Next, in the pure 
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Yang-Mills theory which we consider here, the stochastic vacuum model is only applicable 
up to the distances where the glueball production starts. This happens at |x| 



A/glucball 



^fund 



Setting for an estimate M gme baii = lGeV, we get |x| max = 5.56 GeV 1 . (This yields 
^adjMmax = H-ll, that still leaves lot of space for the inequality o^x 2 3> 1 to hold.) 

(9<r 2 |x|) 1/3 

Accordingly, - — — < 0.39, that demonstrates the numerical smallness of the term 
(9/32) 1 / 3 (cr a djX 2 ) 2 ' /3 with respect to the term a^L\x\. 

Let us finally consider the pre-exponent of eq. (2.4). First, the relation a^L = Mq 
yields L = 0.92 fm. This value for the mean distance, by which a valence gluon is separated 
from an infinitely heavy adjoint source, well exceeds the corresponding vacuum correlation 
length Mq 1 ~ 0.15 fm. The string in a gluelump is therefore a clearly shaped object. One 
can further see that, at this value of L, the term * /3 in the brackets is always larger 

than the term ^ ( "ET ) • Indeed, to compare the latter to the former means to compare 
fsgp = 2 6.08GcV-i ; tQ 2/3 (at |x| » ^ = 1.67 GeV" 1 ). Clearly, 2 -^f^ < f 
up to | a; | = 39.12 GeV -1 , that is the distance much larger than |x| max . This proves our 

1 /3 

statement that 1 ^ L/3 > ^ • As f° r ^ ne term 3^ .^2)2/3 1 it exceeds only in 

the limiting case cr a ^x 2 < 1.89, that is on the border of applicability of the inequality 
Cadj^ 2 S> 1. Thus, we conclude that, in general, the term ^ 1/3 is the defining one for the 
pre-exponential behavior of the function D\(x). This behavior is, therefore, ~ x~ 2 , that 
differs from the behavior ~ |a;| _1 of eq. (E 



3. Green function of a two-gluon gluelump and the function D 



This gluelump is formed by two gluons, which are attached by fundamental strings to an 
infinitely heavy adjoint source and joined to each other by a third fundamental string. 
Since the gluons are identical, the vectors of their spatial coordinates, x(r) and y(r), can 
be chosen in the form x(r) = L'ni+li(r), y(r) = L'xi2+h{T)- Here, ni = (cos 9, — s'mO, 0), 
112 = (cos 0, sin 0,0), and 8 is some constant angle, whose possible values as a function of 
V will be discussed at the end of the calculation. The mean distance L' from the gluons 
to the heavy source is, in general, different from L. Within such notations, the area swept 
out by the string joining the gluons is 2L'T sin9. After the extraction of the mean-area 
dependent factor, the scalar function G^ 2gl \T, 0), through which the Green function of the 



gluelump is expressed [10|, factorizes into the product of integrals over trajectories of each 
of the gluons: 



G (2gl) (T, 0) 



-2(Tf UIK jL'T sin t 



ds 




ds' J VX4 

x 4 (0)=0,x 4 (s)=T y 4 (0)=0,s/ 4 (s)= 



V Vi / PxPyx 



x exp 



x 2 

dr y— + -1 + T] U „,i|x|.r 1 



y2 y2 

dr ( — + — + <7 fund |y|y 4 



(3.1) 
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Let us consider, for example, the integral over x(r). Substituting |x| ~ V + nilj, we have 
for it 

4 2 2/x 

We can further use the same approximation as in the previous section, namely substitute 
(nili) 2 — > ^ J"q c?r(nili) 2 . Together with the (free) integration over lf(r), this yields 



j V\\D[i exp 


KM 







\747rs 



VTZV[i exp 



where 7£ = (^,^ 2 ) and 

E/Yft) = f ^Li 1 {l\ cos e-l\ smO f. 

Using the variational method in the 2-d case (see p| for details), we further obtain 
j VKexp 



4 



1 

Airs 



U{K) 



exp 



d 2 KU(K)e- 



(7 



fund- 

12 



r 2 



The resulting integral over x^ and /i then coincides (modulo the replacement o" ac jj — > ffund) 
with that of eq. (p.3|) , so that the total integral over efePa^Dx in eq. (3.1) reads 



(afund/T 4 ) 1 / 3 
3 2/3 . 2 7/3 . ^3/2 



exp 



-a fund L'T-(9/32) 1 /3( c7fund r 



^2\2/3 



The integral over ds'T>y/{Dy apparently produces the same result. Taking into account the 
relation @ D(x) = ^f^G (2gl) (x, 0), we finally obtain that, at cr fund £ 2 > 1, 



D{x) 



g\N 2 c 



( 0"fund 

I - — — I exp 



3 4/3 . 2 17/3 . ^3 ^ \ x \4 J 



-2a fund L'(l + sin0)|x| - (9/4) 1 / 3 (a fund x 2 ) 2 / 3 . 

(3.2) 

Again, as in the case of the function D±(x), the term 2<7f un d£'(l + sin 0)\x\ in the exponent 
is numerically larger than (9/4) 1 / 3 (a"f un dX 2 ) 2/ ' 3 . Indeed, identifying 2<7f un d£'(l + sin 6) with 



M^ 2gl) = 2.56 GeV from the fori mil;. |1() 



D{x) „ 9 (^1) 0-1Q8(T 2 ^e-J^N va ii d at M ( 2gl) |x| » 1, 



; fund 



(3.3) 



WP have (9/ 4 ) 1/3 (^n d x 2 ) 2 / 3 < (9< nd NUax/4) 1/3 
We haVe 2 CTfund L'(l+sin 6)\x\ ^ ^p) 



0.12. Furthermore, the mean distance 



between the two gluons, as well as the mean distance between these gluons and the heavy 



U < \x\ 



should 



source may not exceed | x| max , i.e. the inequalities 2.Z7sin$ < l^-lmax 
be respected together with the equation 2crf un d^'(l + sin#) = M^ gl \ One can readily see 
that the unphysical case 6 = cannot be realized also numerically. Indeed, one should 
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then have — = 2L' < 2\ x| maX) that is not feasible, since the l.h.s. equals 14.22 GeV , 

CT fund 

whereas the r.h.s. equals 11.12 GeV -1 . For all other cases, we have: 9 < arcsin ^ 
at Mpi < J] < | x | max , while no restrictions on the values of 9 exist at L' < ^^ iax . 3 
Finally, we note that the pre-exponent in the function fl3.2|) falls off as |x| -8 / 3 , whereas the 
pre-exponent of the function (|3.3| ) is a constant. 



4. Summary 

In the present paper, we have evaluated the functions D(x) and D\(x), which parametrize 
the two-point gauge-invariant correlation function of gluonic field strengths in the stochastic 
vacuum model of QCD. To this end, we have applied the Feynman variational method to 
the calculation of the Green functions of one- and two-gluon gluelumps, through which the 
functions D and D\ can be expressed. At the end of the calculation, we have compared our 
results with those found by the Hamiltonian analysis in ref. [pi]. In case of the function 
D±, this comparison yielded the mean separation of a gluon from the heavy source in the 
one-gluon gluelump. In case of the function D, it yielded some constraints on the values 
of the angle between the planar world sheets of the gluons in the two-gluon gluelump. We 
have further seen that, although fluctuations of the gluons (which are accounted for by the 
variational method) around their mean separations from the source produce the terms in 
the exponential asymptotics of the functions D and Di, which are parametrically larger 
than the leading ones, numerically these terms are nevertheless subleading, as they should 
be. Furthermore, we have also found the distance dependences of the pre-exponential 
factors in the both functions D and D%, which turn out to be different from those of 
ref. JlOfl . These differences are, however, of a rather academic interest, since numerically 
a distance dependence in the pre-exponent plays an unimportant role at large distances 



considered both here and in ref. [1C] 



In conclusion, the obtained results might be important for a field-theoretic establish- 
ment of the stochastic vacuum model. Besides that, the performed calculation demon- 
strates how one can account for the fact that the virtual particles in QCD diagrams are 
confined rather than free. Last but not least, the novel expressions for D and D\ obtained 
here might be useful for applications in high-energy scattering. 
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